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Abstract 

We consider two-species exclusion processes on the d-dimensional discrete torus taking 
the effects of exchange, creation and annihilation into account. The model is, in 
general, of nongradient type. We prove that the (charged) particle density converges 
to the solution of a certain nonlinear diffusion equation under the diffusive rescaling 
in space and time. We also prove a lower bound on the spectral gap for the generator 
of the process confined in a finite volume. D | 

Keywords: hydrodynamic limit; interacting particle systems; two-species exclusion pro- 
cesses 



1 Introduction 

The aim of this paper is to obtain the hydrodynamic behavior of two-species exclusion 
processes. Our results can be applied to establish the hydrodynamic limit for the evolution 
of height differences in interfaces governed by the 1-dimensional SOS dynamics. 

The two-species exclusion process describes the evolution of a system of mechanically 
distinguishable particles, say -|-particles and —particles moving on a discrete lattice space 
under the constraint that at most one particle can occupy each site. The state space of 
the process is given by {—1,0, 1}'^n where stands for the d-dimensional discrete torus 
with side-length N and its elements (called configurations) are denoted by 77 = {r]{x),x £ 
T^), with r/(x) = or 1 or —1 depending on whether x G is empty or occupied 
by a -l-particle or a —particle, respectively. Each it particle moves to a neighboring 
empty site with the constant jump rate C± > 0, respectively. Two different types of 
neighboring particles exchange their locations with the constant rate C^; > 0. Also they 
annihilate simultaneously when they are neighboring with the constant rate Ca > 0, and 
two different types of particles are created with the constant rate Cc > if two empty 
sites are neighboring. 

In this paper, we consider the case where Ca > or Cc > 0, so that the process has a 
unique conserved quantity "^^^jn rj{x). We prove the hydrodynamic limit for the profile 
associated with this quantity and obtain the explicit expression of the diffusion coefficient. 
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We classify the dynamics into three types as the case where Ca > and Cc > (Case 
1), Ca > and Cc = (Case 2) and Ca = and Cc > (Case 3) and give proofs 
separately. For Cases 2 and 3, we assume the gradient condition so far. We can show that 
all of the hydrodynamic equations of our processes have a diagonal diffusion coefficient 
matrix, therefore they have unique weak solutions without the smoothness of the diffusion 
coefficient. 

Quastel proved the hydrodynamic limit for two-colored simple exclusion process in [6]. 
This process is corresponding to the model with C+ = C_ > and Ca = Cc = Ce = 0, 
so it is not included in our cases. 

The SOS dynamics describe the evolution of the integer-valued heights of interfaces 
on the discrete lattice. In the 1-dimensional case, the height difference of SOS dynamics 
and the configuration of the two-species exclusion process have one-to-one correspondence, 
see, e.g. [J. 

This paper is organized as follows: In Section 2 we introduce our model and state the 
main results for three types of models respectively. In Section 3, we give the proof of the 
main theorem for Case 1. In the proof, we give a spectral gap estimate and characterize 
the class of closed forms. In Sections 4 and 5, we give the proofs of the main theorems 
for Cases 2 and 3, respectively. In Section 6, we state the uniqueness results for nonlinear 
parabolic equations whose diffusion coefficient matrices are diagonal. 

2 Model and Main Results 

The two-species exclusion process is a Markov process r]t on the configuration space x'n = 
{—1,0, l^^N , where = (Z/A^Z)'^ is the d-dimensional discrete torus. The dynamics are 
defined by means of an infinitesimal generator Ljq acting on functions / : x% ^ K as 

(L^/)(r/)= ^^■/■(^)' 

where (T^)* stands for the set of all directed bonds b = {x,y), i.e., the ordered pairs of 
x, y S such that \x — y\ = 1 where |x — y| = '}2i<i<d ~ y«l sxvm norm in M'^. 

Here, for each bond h G (T^)*, 

(2.1) W(r/)=Cb(r/)(7rb/)(r/), 
where 

i'^bDiv) = [i{rj(x)=l,r,{y)=0} + l{»7W=0,»?(j/)=-l} + l{»)(x)=-l,r;(y)=l}] (/(^^'^) " /(??)) 

+ l{,(,)=i,,(,)=_i}(/(r/-'=°'^=°) - /(r?)) + l{,(.)=o,,(,)=o}(/(^"=~'''^=') - fm, 

and 

Cfe(??) = C+l{^(^.)=i^^(y)=o} + C'_l{^(a,)=o,,y{y)=-l} + CEl{rj(x)=-l,Ti{y)=l} 
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In the above formula, rf'^^rf ^'^ ^ and rf ^ G x% stand for 



r]{z) \iz^x,y, 
v{y) if z = x, 
^r]{x) iiz = y, 



f=-^^y=\z) = 



r]{z) iizy^x,y, 
— 1 ii z = X, 
1 a z = y, 



and 



n 



ri{z) iizj^x,y, 
ii z = X, 
iiz = y, 

respectively, and we assume that C+ and C_ are positive constants and Ca,Cc and Ce 

are nonnegativc constants. 

We will use the following simplified notations 



7]- 



x.y 



if 77(y)) = (1, 0) or (0, -1) or (-1, 1), 



?7"=°'^=° if (??(:r),r?(y)) = (1,-1), 
y^x=-l,y=l if (^(a;),,,(y)) = (0,0), 



and define for 6 G (1"^)* by 

With these notations, Lf, and tt^ can be rewritten as 

W(^)=C5(r/)(/(r/(^'^))-/(7?)), 
^6/(^) = r6(7?)(/(77("'^')) -/(??)), 

respectively. 

The process is reversible with respect to the following one parameter family of trans- 
lation invariant product measures fp. 

Definition 2.1. For each fixed p G [—1,1], let Up be a product measure on Xn 
marginals given by 



i^pivix) = 1} 



1 - + p 



i^pMx) = 0} = 



Mv{x) = -1} 



1 - - p 
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for all X G T^, where 



Hp) = r '-'^ 



with 13 = Cc/Ca- Especially, if Ca > 0,Cc = 0, then 

Up{7]{x) = 1} =pVO, 
UpMx) = 0} = I - \p\, 
i^,{r]{x) = -1} = \pAO\, 



and if Ca = 0, Cc > 0, then 





= 1} 




= 0} 




-1} 



l±p 

2 ' 

0, 

l-p 



The index p stands for the density of particles with charge, namely Ei^p[ri(0)] = p. We 
will abuse the same notation Up for the product measures on the configuration spaces Xn 
or x'^ = {~1)0, 1}^ on the torus or on the infinite lattice. The expectation with respect 
to Vp will be sometimes denoted by 

' f{rj)up{dT^) = {f)p. 

Prom the definition, our model satisfies the detailed balance condition, namely, for 
any directed bond b = 

(2.2) c,(r?)z.,(r?)=c,K??^"'^Vp(^^"''^) 

holds, where b' = (y, x) is the reversed bond of b. 

Here and after, we call f a cylinder function on x'^ if f depends on the configurations 
only through a finite set of coordinates. For any directed bond b = (x, y) and cylinder 
functions /, g, let us define 'Db{i>p; f, g) and Vi,{vp;f) by 

T^b{vp;f,g) ■■= {-iLb + Lb')f,g)p 

and 



where b' = {y,x) and (•, ■)p stands for the inner product in L^(fp). The reversibility ()2.2 



Vb{up;f) ■.= Vb{up;fJ), 

n L'^ihpj 
implies 

(2.3) Vb{up;f,g) = {cb{-Kbf){T^bg)) p- 

Let Tx be the shift operator acting on the set A dll^ and cylinder functions / as well 
as configurations r/ as follows: 
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TxA:=x + A, Txf{ri) = f{Txri), {Txr]){z) := r]{z - x), zeZ'^. 
For every cylinder function g : consider the formal sum 



which does not make sense but for which the gradient 

T^Tg = (7ro,eir(„...,7ro,erfr<;) 

is well defined. 

We are now in a position to define the diffusion coefficient. For each p G [—1,1], 
define 

d{p) = ^ inf Po,e(z^p; m + Tg) 

where inf^ is taken over all cylinder functions g and e is a unit vector of arbitrary direction. 
In this formula x(p) stands for the so-called static compressibility which in our case is equal 
to 

x{p) = m')p - {riml = 1 - '^(p) ~ p' 

Notice that d{p) does not depend on the choice of a unit vector e^, 1 < i < d. 

For a probability measure /x^ on Xad we denote by P^jv the distribution on the path 
space -D(M+, Xn) of the Markov process rjt = {r]t{x), x G T^} with generator N'^Ln, which 
is accelerated by a factor iV^, and the initial measure /x^. Hereafter E^jv stands for the 
expectation with respect to P^iv. 

With these notations our main theorems are stated as follows: 

Theorem 2.1. Assume Ca > and Cc > 0. Let {ij-^)n>i be a sequence of probability 
measures on xff such that the corresponding initial density fields satisfy 

lim /.^[l-^ J2 Gi^Mx) - [ Giu)po{u)du\ > S] = 0, 

for every 5 > 0, every continuous function G : T*^ — M and some measurable function 
po :T'^ [—1, 1]. Then, for every t > 0, 

limsupP^iv[|-^ J2 G{^)Vt{x) - [ G{u)p{t,u)du\ > S] = 0, 



for every 6 > and every continuous function G : T'^ ^ M., where p(t, u) is the unique 
weak solution of the following nonlinear parabolic equation: 



(2.4) 



where 



dtp(t,u) = A(J(p(t,«)))( = ^£-{d(p(i,«))^(i,«)}) 

i=l 

p{0,-)=Po{-), 



d{p) = j\^)d^. 
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The rigorous definition of weak solutions is given in Section 6. 

Remark 2.1. // we assume C+ + C_ — Ca — "^Ce = 0, then our model turns out to be a 
gradient system. In this case, d{p) = - *M(c^ _ c_) + 1 (c+ + C_) holds. In particular, 
we can compute the diffusion coefficient d{p) explicitly from the concrete values ofC^,C- 
and (3. 

Remark 2.2. Generalized exclusion process with k = 2 is corresponding to our model with 
C+ = C- = Ca = Cc = 1 and Ce = 0. 

Theorem 2.2. Assume Ca > 0, Cc = and the gradient condition C++C- — Ca — 2Ce = 
0. Let {fJ.^)N>i satisfy the same assumption as in Theorem \2.1l Then, for every t > 0, 

limsupP^iv[|-^ V G{^)r]tix) - [ G{u)p{t,u)du\ > 5] = 0, 

for every 5 > and every continuous function G : T'^ ^ M, where p{t, u) is the unique 
weak solution of the following nonlinear parabolic equation: 




{ p(0,-) =/5o(-), 



where 

(2.6) P{p) = C+pl|,>o} - C_plv<o}- 

Remark 2.3. Equation ^2. 5\) is the weak or enthalpy formulation of the following two- 
phases Stefan problem: 

' dtp{t, u) = C+Ap{t, u) on C{t) = {p{t, u) > 0} 
dtp{t, u) = C-Ap{t, u) on S{t) = {p{t, u) < 0} 

< = n • (^C+V{p{t, u) V 0) - C_ V(p(t, u) A 0)) 

on T,{t) = {p{t,u) = 0} 

where n denotes the unit normal vector on directed to C{t) and V(/>(t,n) V 0) (re- 
spectively V{p(t,u) AO)) is the limit of the gradient o/pVO (respectively pAO) at u £ T,{t) 
when approached from C{t) (respectively S(t)), see 12]. 

Theorem 2.3. Assume Ca = 0, Cc > and the gradient condition C+ + C_ — 2Ce = 0. 
Let {p^)n>i satisfy the same assumption as in Theorem \2.1\ Then, for every t>0, 

limsupP^iv[|-^ V G(^)r/t(x) - / G{u)p{t,u)du\ > 5] = 0, 



6 



for every 5 > and every continuous function G : T"^ ^ M, where p{t, u) is the unique 
weak solution of the heat equation: 



(2.7) 



dtpit, u) = CEApit, u)( = CeY, fT^P^^^ 



■ 1 

1=1 ' 



P(0,-) =/Oo(-)- 
3 Proof of Theorem 12.1 



In this section, we consider Case 1, namely the dynamics with Ca > and Cc > 0. 
For this case, we do not assume the gradient condition, so the system is, in general, 
nongradient. The strategy of the proof is essentially the same as given for the generalized 
exclusion process in The main step is obtaining the estimate of the spectral gap 
and the characterization of the closed forms, which are presented in the Subsections 3.4 
and 3.5. To guarantee the uniqueness of the weak solution, we show that the diffusion 
coefficient matrix is diagonal. The method of the proof is also available for the large class 
of symmetric processes including generalized exclusion process. 

3.1 The Macroscopic Equation 

We start with considering a class of martingales associated with the empirical measure. 
We take T > arbitrarily and fix it in the rest of this section. For each smooth function 
H ■.T'^ ^R, let M^'^{t) = M"{t) be the martingale defined by 

M^(t) = (^f - (vTo^,/?) - f N'LN{7Tf,H)ds, 

Jo 



where nj- stands for the empirical measure associated with rjt, namely 

(3.1) vrf (dn) = ^Y1 Vtix)S^idu), 0<t<T, uGT'', 

and (vr^, /) stands for the integration of / with respect to vr^. 

A simple computation shows that the expected value of the quadratic variation of 
M^[t) vanishes as | oo, and therefore by Doob's inequality, for every 5 > 0, we have 

lim P iv[ sup \M"{t)\ > (5] = 0. 

N^oo 0<t<T 

A spatial summation by parts permits to rewrite the martingale M^{t) as 
M^{t) = {nf,H)-{n^,H) 

(3.2) - X] fN'-" Yl {dZH){^)r.W,^eAris)ds, 
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where Wo.Ci represents the instantaneous current from to e^: 

+ C_(<l\(r?)-^°j!:o(^)) 
and d^.H represents the discrete derivative of H in the i-th direction: 

Next we show that the current Wo,e- can be decomposed into a hnear combination of 
the gradients {rj^ej) — ??(0), 1 < i < d} and a function in the range of the generator L^v: 
^o,ei + Z]j=i ~ "^(0)] = -^Aff foi' ^ certain cyhnder function f and some matrix 

Dij(p) that depends on the density, see Theorem 13.11 and Corollarv 13.21 for more precise 
statement. 

Denote by {Di j(p), 1 < i, j < d} the unique symmetric matrix such that 

1 

a*D{p)a = ^^inf { V2?o,e,(i^p, ai??(0) + T,,) |, 

for every vector a in where inf^ is taken over all cyhnder functions g. 

For positive integers I, N, a function H in C^(T'^) and a cylinder function f on x'^i let 



where 



and 



(r?) = Tyo,e,(r?) A,,(r/^(0))[r/'(e,) -r?'(0)] -LArf(r?), 



iV- 



Theorem 3.1. Fix p G (—1,1) arbitrarily. Then, for every function H in C^{T'^) and 
1 < i < d, we have 

1 

inflimsuplimsup— -TlogE^iv[exp{A^'^| / x]^ j^{H,r]s)ds\}] = 0, 
fee Af^oo -'v" Jo ' 

where C stands for the set of cylinder functions on x'^- 

The proof of Theorem 13. II is postponed to the next subsection. This theorem imphes 
the following corollary. For a positive integer / and a function H in C^(T'^), let 

d 
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Corollary 3.2. For every function H in C^(T'^) and \ <i < d, 

r ■ 

limsuplimsupE^]v[| / Y}^^^^{H,r]s)ds\] = 
e^O N-*oo Jo 

To prove this corollary, we can use the method in [3] straightforwardly. In particular, 
the L]\ff term is negligible. We have now all elements to prove the hydrodynamic behavior 
of our nongradient system. 



Proof of Theorem \2.1\ Recall that the empirical measure vr^ is defined by (|3.1|) . Denote 



by Q^N the distribution on the path space £)([0,T], A^(T'^)) of the process vr^ where 
M{T'^) stands for the space of signed measures on T'^ endowed with the weak topology. 

Following the same argument as for the generalized exclusion process in [3] it is easy 
to prove that the sequence {Q^NjN > 1} is weakly relatively compact and that every 
limit points Q* is concentrated on absolutely continuous paths irtidu) = Tr{t,u)du with 
density bounded by 1 and -1 from above and below respectively: —1 < TT{t,u) < 1. 

From Theorem 16.11 stated below, there exists at most one weak solution of (j2.4p . 
Therefore, to conclude the proof of the theorem, it remains to show that all limit points 
of the sequence {Q^n,N > 1} are concentrated on absolutely continuous trajectories 
7r{t,du) = 7r{t,u)du whose densities are weak solutions of the equation (j2.4|) . 

Fix a smooth function : T'^ ^ M and recall the definition of the martingale M^{t). 
Applying Corollary 13.21 to the last integral term in the formula ()3.2p of M^{t), we obtain 
that for every 6 > 0, 

limsuplimsupP ]v[|(7r^,i7) - {it^,H) 

where 

Denote by Dij the integral of Dij : Dij{p) = J^^ Dij{'~f)d'y. Since H is smooth and Dij 
is continuous by Theorem 13.81 stated below, with the help of Taylor's expansion and a 
spatial summation by parts, we have 

limsuplimsupP iv[|(7r^,iJ) - {Tr^,H) 

s^O N^oo 

Therefore, for every limit point Q* of the sequence Q^n , 
limsupg*[|(7rT,ii')-(^0,^) 



f ds j duidl^^M){u)Di,,{{7rs * t,){u))\ > 6] = 
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where 

:= (2e)-^lf_,,,].(.) 

and * represents the convolution. Since each hmit point Q* is concentrated on absolutely 
continuous paths nt = vr(t, u)du with —1 < 7r(t, u) < 1, for each fixed < s < T, {'Ks*i^e){u) 
converges to 7r(s,u) for almost u in T*^ as e | 0. From this remark and the continuity of 
{Dij, 1 < i,j < d}, we obtain that 



Q*[ {7rT,H) - {7ro,H) - V ds [ du{dl H){u)Diji 



■n{s,u)) 



>5\ 



for all H in C^(T'^). The fact D{p) = d(p)I, namely D{p) = d{p)I proved in Theorem 
13. 101 permits to rewrite the last expression as 

Q*[{7rT,H) = {Tro,H) + [ ds [ duAH{u)diTT{s,u))] = 1. 

Jo Jt^ 

Denote by {tnjngN a dense subset of [0, T] and repeat the same argument as we have done 
up to this point for any fixed t^, then 

Q*[{TTt„, H) = {ttq, H) +[" ds [ duAH{u)d{7:{s,u)) for every n G N] = 1. 

Jo Jjd 

Since Q* is the probability measure on D space and d is bounded function on [—1, 1] and 
Q* is concentrated on paths vrj = Tr{t,u)du with —1 < 7r(t, n) < 1, Q* is concentrated on 
the weak solution of (12. 4p which concludes the proof of the theorem. □ 



3.2 Central Limit Theorem Variances 



To state the main theorem of this subsection, first we introduce some notation. For a fixed 
positive integer / we denote by A; a cube in Z"^ of side-length 21 + 1 centered at the origin: 
A; := {— /, — / + 1, — 1,1}'^. We denote the set of cylinder functions on x'^ by C. For 
^ in C, denote by A^ the smallest d-dimensional rectangle that contains the support of 
and by sq, the smallest positive integer s such that Aij, C A^. Let Cq be the space of 
cylinder functions with mean zero with respect to all canonical invariant measures: 

Co = { 5 e C ; {g)A^,K = for all - jA^I <K< \Ag\ }. 

Here, for a finite subset A of Z*^, we denote by |A| the cardinality of A and by {■)a,k the 
expectation with respect to the canonical measure iyA,K '■= ■ IZ]xgA^(^) ~ ^) 
— |A| < K < |A| which is indeed independent of the choice of a. For a rectangle A and 
a canonical measure va,k, denote by {■,-)a,k (i'esp.(-, the inner product in L?'{va,k) 
(resp. L'^iva))- 

It is known that to conclude the proof of Theorem 13.11 it is enough to show that 
(3.3) inf hm sup(2/)'^((-LAj-iy/'', T//'');,^ = 
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where 

y/''(r/) = (2r + ir'^ j;r,T^O,e,(r/) 

\y\<i' 

+ ^ A,,(r?'(0))[/(e,) - /(O)] - (2^f + 1)"'^ (^v^Nm, 

I' = I — 1 and = / — — 1 so that TyL^^ is ^A;-™6asurable for every y in A^^. This 
follows from Theorem 13.41 and Corollary 13.91 below. 

For the beginning of the proof we obtain a variational formula for this variance. We 
start with introducing a semi- norm on Cq, which is closely related to the central limit 
theorem variance. For 1 < k < d denote U'' = {U^ , the d-dimensional cylinder 

function with coordinates defined by 

iU^Uv) = 6i,k^o,e,m for all I < i < d. 

Here 5ij stands for the delta of Kronecker. For cylinder functions g, h in Cq and 1 < i < d, 
let 

9,h:>pfi= 'Y{g,T^h)p and <^ g :^pj= Xj{g,r]{x))p, 

where xj stands for the j-th coordinate of x G Z'^. Both <^ g,h ^p,o and <^ g are 
well defined because g and h belong to Cq and therefore all but a finite number of terms 

vanish. For h in Cq, define the semi- norm <^ h^p by 
(3.4) 

d d d 

= sup {2<5,/i>p,o+2j]]ai</i>p,, -^((^aj(^/-')i + Vo,e,rg)^)p} 

d d 

= sup {2<C5,/i>p,o+2^ai</i>p,i-^Po,e,(i^p;ai??(0)+rg)}, 
where a = {ai)f^i. 

1 

We investigate in the next section several properties of the semi-norm <^ ■ ^p , while 
in this section we prove that the variance 

\x\<l^ \x\<l^ 

of any cylinder function ^p in Cq converges to <^ ip as / j oo and p. Here 

stands for / — so that the support of T^ip is included in for every x < l^. By 
elementary computations relying on an adequate change of variables, the norm <C • ^p 
may be rewritten as 



d 

^h^p= sup {2^g,h >p,o +2 ^ ai < /i >p,, 

gGCo,aGK'* j=i 
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d 

+ 2 J^a, « Tyo,e„9 »p,o ((Vo,eir/(0))2), - {\\VTgf)p}. 

i=l 

We are now in a position to state the main result of this section. 

Proposition 3.3. Consider a cylinder function ip in Cq and a sequence of integers Ki 
such that -{21 + l)'^ < Ki < {21 + l)'^ and hm^^oo (fp = P- Then, 

/ — 5-00 ^ ' ^ ' 

\x\<l^ \x\<l^ 

Once Theorem 13. 181 which is stated below, is established the proof of Proposition 13.31 
is the same as that of Theorem 7.4.1 of ^ since the proof does not depend on the specific 
form of Pfc. 

We conclude this section proving that for each ij) in Co the function <C V' 1] ~^ 

M_i_ that associates to each density p the value <C ip is continuous and that the 
convergence of the finite volume variances to <C • is uniform on [—1,1]. For each / in N 
and-(2/+l)'^ <K < (2/+l)'^, denote by V;'^((2jfjjd) the variance of (2/+l)~"' ^l^,|<^^ r^V' 
with respect to ui^k- 

^ ' \x\<l^ \x\<l^ 

We may interpolate linearly to extend the definition of V^^ to the all interval [—1, 1]. 
With this definition V^^ is continuous. Proposition [3?3l asserts that Vj^ converges, as / j oo, 
to ^ V for any sequence Ki such that ^2;+i)d ~^ P- particular, lim^^oo Vi^{pi) =<C 

ip for any sequence pi — > p. This implies that <C V' is continuous and that V^'^(-) 
converges uniformly to <C as / | oo. We have thus proved the following theorem. 

Theorem 3.4. For each fixed h in Cq, <^ h is continuous as a function of the density 
p on [—1,1]. Moreover, the variance 

\x\<lh \x\<lh 

converges uniformly to <^ h >>p as / t and ~^ ^' particular, 

lim sup (2/)~'^((-LaJ"^ r^h, T^h)i^K = sup < >p . 

'^'^-i2i+ir<K<{2i+ir -i<p<i 

3.3 The Diffusion Coefficient 

We investigate here the main properties of the semi norm ^ • S>p introduced in the previ- 
ous section. We first define from <^ ■ ^p a semi-inner product on Cq through polarization: 

(3.5) <5,^>p= ^{<5 + ^>p- <5-^>p}- 
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It is easy to check that (I3.5p defines a semi-inner product on Cq. Denote by Mp the 

kernel of the semi-norm <C • 3>p on Cq. Since <C • is a semi-inner product on Cq, the 
completion of ColATp, denoted by Hp, is a Hilbert space. 

Simple computations show that the linear space generated by the currents {VFo.e, ) 1 < 
i < d} and LCq = {Lg; g £ Cq} are subsets of Cq. The first main result of this section 
consists in showing that Tip is the completion of LCqIj^^ + {Wo,ei,l < ^ < d}, in other 
words, that all elements of Tip can be approximated by X]i<i<d ^i^o.cj + Lg for some a 
in M.'^ and g in Cq. To prove this result we derive two elementary identities: 

(3.6) <. h, Lg >p= - <.h,g >p,o and < h, Wo^e, >p= - < /i >p,i 

for all /i, 5 in Cq and 1 < i < d. 

By Proposition 13.31 and (|3.5p . the semi- inner product <^ h,g ^p is the limit of the 
covariance (2/)~'^((-LaJ"^ S|x|</g "^i^ff' Z]|x|<«h '^xh)i,Ki as / j oo and p. In partic- 

ular, if 5 = i(70) for some cylinder function go, the inverse of the generator cancels with 
the generator. Therefore, <C h, Lgo ^p is equal to 

- lim(2/)"'^( V T^go, V Txh)i^Ki 9o, h :^pfi . 

I — >oo — ' ^ — ' 

The second identity is proved in a similar way. 

It follows from the first identity of (13. 6|) that the gradients {ri{ei) — r]{0), 1 < i < d} are 
orthogonal to the space LCq, while the second identity permits to compute inner product 
of cylinder functions with the current: 

(3.7) «r?(e,)-?/(0),L/i»p = 0, 

(3.8) < r]{ei) - r/(0), Wo,e, >p = -xip)kj^ 

and 

(3.9) < T^O,e.,W^O,e, »p= ((Vo,eir?(0))2)/,,,-. 

for all 1 < i, j < d and /i G Cq. In this formula x{p) stands for the static compressibility 
and is equal to (??(0)^)p — (??(0))p. Furthermore, 

d d 

(3.10) < Xl'^i^O'^. + ^9 »p= X]^{^o,e,(air/(0) + Tg)}^)p 

j=i i=i 

for a in M*^ and g in Cq. In particular, the variational formula for <^ h S>p writes 

d d 

(3.11) </i>p= sup {-2<^h,J2aiWo,e,+Lg:^p- 't:^a,Wo,e,+Lg:^p}. 

Proposition 3.5. Recall that we denote by LCq the space {Lg; g £ Cq}. Then, for each 
— 1 < p < 1, we have 

Hp = LC^\x, © {Wo,e,,l <i<d}. 
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Proof. We can apply the proof of Proposition 7.5.2 in [3j straightforwardly. □ 

Corollary 3.6. For each g G Co, there exists a unique vector a G R'^ such that 

d 

5 - ^ ajWo^ej G LCo in Tip. 
i=i 

We now start to describe the diffusion coefficient D of the hydrodynamic equation. 
From Corollary 13. 6^ there exists a matrix {Qij, ^ < i, j < d} such that 

d 

(3.12) T]{e,) - r/(0) + J2 Q*,j^o,e, eLC^ in Hp. 

i=i 

Notice that the matrix Q = Q{p) depends on the density p because the inner product 
depends on p. It is easily shown that Q is symmetric and strictly positive. 

Denote by D = D(p) the inverse of Q, which is also symmetric and strictly positive. 
We will see below that D{p) is the diffusion coefficient of the hydrodynamic equation (|2.4|) . 
Since D is the inverse of Q, we have that 

d 

Wo,e, + J2 DiAviej) - vm G in Hp. 

i=i 

for 1 < i < d. This relation provides a variational characterization of the diffusion 
coefficient D. Indeed, for all vectors a G M*^, 

d d 

(3.13) inf {« V a^Wo,e, + V aiAj[7?(ej) - 7?(0)] - Lg >p} = 0. 

1=1 *)J = 1 

Since gradients are orthogonal to the space LCq, 

< r]{ej) - r/(0), Wo,e, »p= -x{p)k,v 

and 

< ??(ej) - r/(0),r?(efc) - r/(0) »p= x(/5)Qj,fc = x(p)[^~^]i,fc, 
the last identity reduces to 

d 

inf {-x(p)a*Da+ « VaiVFo.ei - ^5 »p} = 0, 

9GC0 

where a* stands for the transposition of a. We have thus obtained a variational formula 
for B{p). 

Theorem 3.7. The diffusion coefficient D{p) is such that 

^ d 1 

(3.14) a* Da = ^ mf^ « ^ a,Wo,e, - Lg »,= Jnf J^({(Vo,e,(a.r/(0) - r,)}^)^ 
for all a G M'^. 
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The second identity follows from equation (I3.10p . Moreover, this formula determines 
the matrix D since D is symmetric. 

It is now easy to prove the diffusion coefficient is continuous including at the boundary 
of [—1, 1]. From the explicit formulas for x{p)i ((^o,ei^(0))^)p and (^'oei)p) have that 
D(p) converges to C+I as p | 1 and C_/ as p | —1. 

Theorem 3.8. The diffusion coefficient D{p) is continuous on [—1,1]. Moreover it con- 
verges to C+I as p ] 1 and C-I as p [ —1. 

From the continuity of the diffusion coefficient we have 
Corollary 3.9. Let D be the matrix defined in Theorem \3.7\ Then, for each 1 <i < d, 

d 

inf sup < Wo,e, + Y DtApM^j) - rim " ^f(^) >p= 0. 
feCo _i<„<i ^ 

This result together with (j3.3p . the definition of VJ^'' and Theorem 13.41 concludes the 
proof of Theorem 13.11 

We conclude this section proving that the diffusion coefficient D is a diagonal matrix 
and it has the same diagonal component, therefore D{p) = d{p)L 

Theorem 3.10. There exists a continuous function d{p) on [—1, 1] such that D{p) = d{p)I 
and 

Proof. Because of the symmetry of the dynamics, it is obvious that D has the same 
diagonal component. It remains to show that D is a diagonal matrix. 

According to [7], the diffusion coefficient matrix defined by the variational formula 
(|3.14p coincides with the diffusion coefficient matrix defined by the Green-Kubo formula 
based on the current-current correlation function: 

1 r 1 1 

a*D{p)a := | ^ a2((Vo,e.r?(0))2), - - / ^ E^^iWac'^'rMdt} 

where Wa '■= Yli Oil^o,ei- Therefore, we have only to prove that 



poo 



for all i 7^ j. In [S], Kipnis and Varadhan proved some equivalent relation about the 
central limit theorem variance. We can use one of them. It holds that 

roo 

/ [^0,e.e^V,iyo,e,]dt = lim ^ E,^ [Wo,e,T,(/,] 

/n — — — 

where g''^ is a solution of the resolvent equation Xg-'^ — Lcf^ = Wo^ej ■ 
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Denote by 9i the reflection operator with respect to along the direction, namely 
for X G Z'', 6iX = (xi, X2, — + Xj+i, x^). We may extend 9i to configurations 
in x'^ and to functions on x'^ naturally: 

{9,r]){x) :=??(M (^^Z)!??) :=/(M- 

Then, for i ^ j, 

Therefore, since is translation invariant and a product measure, 
Since the map x — > {9iX — 2ei) is a bijection, 

X X 

Thus, Y.X Eu, [Wo,e,T,g{] = for all A. □ 

Remark 3.1. If we assume the gradient condition C+ + C_ — Ca — '^Ce = 0, then 
Wo^e, = Hvi^)) - Hviei)) with h{-l) = C+, h{0) = and h{-l) = -C_. In this case, <. 
Wo,e,,Lg >p= holds. Therefore d{p) = {{V^,eM^)f)p = -^{C+-C-) + \{C++C.). 



3.4 Spectral Gap 

In this section, we prove the spectral gap for the two-species exclusion process on finite 
d-dimensional cubes. For a positive integer N , we denote by iliv the box {1, A^}'^ and by 
the space of configurations {—1,0, 1}^^^ _ ^y^q generator of the two-species 

exclusion process on Vt^q with free boundary conditions: 

x,yenff,\x-y\=l 

where Lxy was defined in ()2.ip . 

For — ir^Afl < K < I^InI, we denote by yN,K the hyperplane {"n'^Ylxenjs, ''l(^) ~ ^} 
and by ij,n,k the product measure on yjy conditioned on the hyperplane yN,K- 

mM-) = ^p{-\ Y = K). 

As in the previous sections, expected values with respect to the measure fJ,N,K are denoted 
by {■)n,k- 



{f)N,K ■■= / f{r])flN,K{dv)- 

In the main theorem of this section we prove that the generator L^^ in IP'{ijli\j^k) has 
a spectral gap of order at least A^^^. 
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Theorem 3.11. There exists a positive constant C , which only depends on the constants 
C+,C-,Ca,Cc and Ce, such that for every positive integer N, every integer —\^n\ < 
K < \^n\ and every function f in L'^{fi]^^K) satisfying {f)N,K = 0, 

{f^)N,K <CN^{-Lnj,f,f)N,K- 

We start with showing that {—Lfij^f,f)iy^K is bounded below by C{—Lfij^f,f)N^K 
with some constant C where Lq^ acting on functions as 

LuNfiv) = Yl ^xyf{ri) 

x,y£Q.M,\x-y\=l 

and 

Lxyfiri) = [^l${ri) + + *"-ti(r/)](/(r?^'^) - f{r,)) 

+ ^?;^i(r/)(/(r/-=°'^=°) - /(t?)) + «:o - fin))- 

Notice that is the generator of two-species exclusion process with C+ = C_ = Ca = 
Ce = 1 and Cc = P- The probability measures fiN,K are also reversible for the Markov 
process with generator Lqj^. 

Lemma 3.12. // we assume that C^,C-,Ca,Cc are all positive constants and Ce is a 
nonnegative constant, there exists a positive constant C such that for every positive integer 
N, every integer — |J7iv| < K < every function f in L'^{^n^k) and every directed 

bond b = {x, y) we have 

{i—Lxy — Lyx)f,f)N,K < C{{—Lxy — Lyx)f, f)N,K 

Proof. It is enough to prove the lemma assuming C^; = 0. Especially, we only have to 
bound the term (^'!.'i^i(r/)(/(r/^'^) - f{v))'^)N,K by the term C{{-Lxy - Lyx)f, f)N,K with 
some constant C. By the Cauchy-Shwartz inequality, we have 

(^^ti(r?)(/(77^'^) - f{v))')N,K = i^llMifiv'''') - f{ri)?)N,K 
< 2{^'l'y_,{v)[{f{v^'y) - /(r/^=°'^=°))2 + (/(r?^=°'^=°) - f{v))'])N,K 

and the last expression is written as 

2/3(^S;o'(^)(/(^'="'''=') - fin))')N,K + 2{^ll,{r^){f{rj^=''y=^) - f{r^))')N,K 

by change of variables. Therefore, we can obtain the desirable estimate with the constant 
C ■.= mm{C+,C-,^}. □ 

Now, to conclude the proof of Theorem 13. IH we have only to prove the theorem as 
follows: 

Theorem 3.13. There exists a positive constant C such that for every positive integer N , 
every — ir^Tvl <K< \^n\ and every function f in LP'{pln,k) satisfying {f)N,K = 0, 

{f^)N,K < CN^{-Ln^f,f)N,K- 
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The proof of this theorem reUes on the study of the spectral gap of the two-species 
exclusion process of mean field type. This is the Markov process on y^, whose generator 
Lq^ acting on functions / as 

Notice that the probability measures hn,k are also reversible for the Markov process with 
generator . This generator has a spectral gap in L^(/xjv,A') of order at least 1 as stated 
in the next theorem. 

Theorem 3.14. There exists a finite constant C such that for every positive integer N , 
every integer < K < \i^N\ o-nd every function f in L'^^^n^k) satisfying {f)N,K = 0, 

{f'^)N,K < C{-L'Q^f,f)N,K- 

Before proving Theorem I3.14t we show that Theorem 13.131 is an easy corollary of this 
result. 

Proof of Theorem \3.13l For each pair {x,y} £ JIa? x O^V! we determine a path inside f^Ar 
which connects x = {xi,...,Xd) and y = {yi,---,yd) as follows: First we connect x and 
{yi,X2, ■■■,Xd) only by changing the first coordinate one by one. Then, (yi, X2, x^..., x^) 
and {yi,y2,X3...,Xd) are connected by changing the second coordinate and this proce- 
dure is continued. We denote the sequence of bonds appearing in this path by bi = 
{zi,wi), b2, bu = {zm,wm) and the set of these bonds by B(x, y). For a configuration t] 
satisfying rx^y{rj) ^ 0, let define a sequence of configurations (Ci)o<j<2A/-i such that = 
C2M-1 = as follows : Co := := (O-i)"^'"^ for j < M - 1, 6/ = {Cm-iY''''""^^ 

and := (^j_;^)^2a/-j,«'2a/-j for M + 1 < j < 2M - 1. Then, we have 

0<i<2M-l 

< (2M - 1) J2 - fiv)?)N,K 

0<i<2M-l 

<MN E {{-Lb-Lb')f,f)N,K. 

beM{x,y) 

Similarly, we have 

{{-Lxy - Lyx)f, f)N,K <CN E {{-Lb - Ly)f, f)N^K 

bm(x,y) 

for some positive constant C . 

Applying Theorem I3.14t for all functions / in L?'{^n^k) satisfying {f)N,K = we 
obtain that 

{f'^)N,K < C—— {-Lxyf,f)N,K<C—-—r {{-Lb - Ly)f, f)N,K 

\\Ln\ — "Af — — 

x,y&\lM x,y€Qi^ b€B{x,y) 
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^^TTTl 1^ {-Lhf,f)N,KX#{ix,y)£nNxnN;borb' eB{x,y)} 
<CN^ ^ {-Lbf,f)A^,K 

feG(njv)* 

where a constant C changes each hne. □ 

Proof of Theorem \3.14\ There is a duahty between +particles and —particles, i.e., +particles 
evolve with the same dynamics as —particles do under the generator Lf^^ . Therefore, we 
assume that < K < \^n\- 

Let X{r]) denote the number of sites occupied by —particles in the configuration rj: 

Xirj) := J2 ^{r,{x)=-i}- 

We first project / on the o"- field generated by X and on its orthogonal: 

(3.15) {f')N,K = ((/ - E[f\X]f)^,K + {{E[f\X]f)^,K. 

We consider the two terms separately. Let us define Ll.y for each ordered pair (x, y) by 

{lyM = i^l'oiv) + ^oUv) + ^^-limiiv''') - f{ri)) 
and define by 

1',{/Gnjv 

To bound the first term in (j3.15p by the Dirichlet form {—Lj^ /, /), we use a general result 
concerning the spectral gap for multispecies exclusion processes. 

We introduce some notation. For positive integers r, and nonnegative integers 
Ki, ...Kr such that Yll=i — define as the hyperplane of all configurations 

of := {0, 1, 2...r}'^ with Ki sites occupied by the z-particles: 

N 

i=i 

and rnj\f,Ki,...Kr as the uniform probability measure on S^^^ As before we denote 

by {') N,Kx,---Kri the expectation with respect to the measure mN^Ki,...Kr- Consider the 
process that exchanges the value of configurations between any two sites at a fixed rate. 
Its generator is given by 

l<j,k<N 

where 

ri(z) if z ^ x,y 
■q^'^{z) = I r}{k) if z = j 

A simple computation shows that the uniform measures raM,Kx,---Kr are reversible for this 
process. We prove that the spectral gap of the generator is of order 0(1). 
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Proposition 3.15. There exists a positive constant C = C{r) such that for every positive 
integer N, every set of nonnegative integers Ki,...Kr such that Yll=i-^i — ^ '^^'^ every 
function f in L'^{m]y^Ki,...Kr) satisfying {f) N,Ki,...Kr = 0, 

(3-16) {f'^)N,Ki,...Kr < C{-LNf,f)N,Ki,...Kr- 

The proof of this proposition is postponed to the last part of this section. To apply 
the estimate in (I3.16P , we rewrite the first term of in the right hand side of ()3.15p as 

l"ivl--ff 

{if-E[f\X]f),,,K= Yl f^NM{X = l}){f?)A^,K,l 

1=0 

where (•)Ajv,-ftr,/ stands for the expectation with respect to the uniform measure on the set of 
configurations rj e {-1, 0, 1}*^^ satisfying Y^xen^ '^Mx)=i} = and Y^x^Qn ^M^)=~^} = 
I, and // stands for the function on this set defined by //(r/) = f{ri) — E[f\X](r]). By Propo- 
sition [3?15l we have that 

l"ivl--g 

{{f-E[f\X]f)N,K < C Yl m,K{{X = l}){-LhjlJl)A^,K,l 

1=0 

= C{-LQj^f,f)N,K 
< C{-L^^f, f)N,K, 

notice that /; can be replaced by / to have the second line. 

Next, let us consider the second term of (j3.15p . Let rjt be the Markov process with 
the generator L"^^. Since the original geometry of the process evolving according to the 
generator L^^ is lost, X{rjt) is a Markov process and the state space of this Markov process 
is X •= {0, 1, ^^^2~^ }- ^ simple computation shows that its generator is given by 

C-N,Kf{l) = r{l, I - !)(/(/ - 1) - /(/)) + r(/, / + !)(/(/ + 1) - /(/)) V/ G x 
where r{l, I - I) = ^^gp and r{l, / + 1) = (\^N\-K-2imM\-K-2l-i)i3 _ ^^^^ ^ 
denote by mN,K the probability measure fiN,KX^^ on x- 

rriN^Kil) ■■= tJ'N,K{{X = I}). 

For X-measurable function /, define a function / : x — > M by /(/) := /(r/) for some r] such 
that X{r]) = I. A simple computation shows that L^f = CN,Kf and {f)N,K = {f)fhjs,K- 

Therefore, {-L%f, f) N,K = {-^N,Kf,f)mN,K and {f'^)N, K = {P )fhN,K hold. To conclude 
the proof of Theorem 13. 141 we have only to prove Lemma 13.161 below. □ 

Lemma 3.16. There exists a constant Cp such that for any integer N and K satisfying 
0<K< \Qn\ 

for all functions / : x ^ 1^ satisfying {f)mMK ~ ^ where rhj^^K, ^n,k o-nd x were defined 
above. 
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The proof of this lemma is based on a general result concerning the spectral gap for 
strongly asymmetric reversible Markov processes presented below, see [1]. 

Proposition 3.17. Let {Xt) be a Markov process with generator denoted by L on a count- 
able state space E reversible with respect to a probability measure m, where L acting on 
functions as Lf{x) = J2yeE^i^^y)ifiy) ~ /(^)) f'^^ x G E. Suppose that there exists a 
point eQ £ E, a positive constant Cq and a ramification {7(60, G E} satisfying the 
following assumption (H): |7(eo,x)| < Co [r{x,p{x)) — Yliyes{x) ^i^iV)]- Then, for every f 
in Lp'{m), we have 

{{f-{f)mf)m<2Co{-Lf,f)m 

Proof of Lemma \3.16l For K = 2\0,i\f\ — 1 and 2|r2Ar| the process has only one possible 
state. Therefore, we may assume that < K < 2\^}j\f\ — 2. 

Denote by the Markov process on x with generator Cn,k- The proof consists in 
finding a state cq and a ramification {7(60, x), a; G Xn} satisfying the assumption {H) 
required in Proposition 13.171 for some strictly positive constant Cq. 

The natural candidate as root of the ramification is the point where the drift of the 
particle is 0. Define D{x) as the mean drift of the particle at x: D{x) = r(x,x + 1) — 
r{x,x — 1). Since < K < 2\0,i\i\ — 2, a simple computation shows that D{0) > and 
D{[ ^^^2 ^^ ]) ^ 0- Let eb be the unique root of D{x) in [0, -^ ]] and define cq as the 
nearest integer of cq. 

Once eo is defined, there is only one possible ramification of the state space. For 
X > Cq we have to define the path from eo to x as 7(eo, x) = (eo, eo + 1, x). In the same 
way for x < eo the path from eo to x has to be 7(eo,x) = (eo,eo — 1, With this 

ramification, for a fixed x > eo the parent of x is x — 1 and there is only one child x + 1. 
In the same way for a fixed x < eo the parent of x is x + 1 and there is only one child 
X — 1. Therefore, in order to prove that this ramification satisfies assumption (H) of the 
Proposition 13.171 we have to show that r(x,x + 1) — r(x,x — 1) > Co(eo — x) for x < eo 
and r(x, x — 1) — r(x, x + 1) > Co(x — eo) for x > eo. 

A simple computation shows that D'{x) < —C^ where = ^^]^^^\^ for (i >\ and 

= 2/? ^^j|^| for 13 <\. Therefore, for x > eo and large N, 

r(x, X — 1) — r(x, X + 1) = —D{x) > (x — e~o) 

> C^{x - eo) - C^leo - e'ol 

>^{x-eo). 

The last inequality follows from the definition of eo since |eo — eo| < ^ and from the fact 
that X > eo + 1. In the same way we can prove that for x < eo, 

r(x, X + 1) — r(x, X — 1) > ^-(eo — x). 

Since there exists some finite constant Cr such that < Cr for all N, we conclude the 
proof of the Lemma 13.161 □ 
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Proof of Proposition \3.15l The idea of the proof consists in using the means of the math- 
ematical induction with respect to r. Assume r = 1. Then, the process is the uniform 
symmetric simple exclusion process for which Quastel proved a spectral gap in [8]. 

Next, consider the general positive integer r. First of all, we suppose without loss 
of generality that Kq < Ki for < i < r where Kq := N — Yll=i^i- This implies that 
Ko < For < i < r, define ttq : {0, 1}^ as a function on the configuration 

space which do not distinguish sites occupied by some particle: 



Mv) = e e {0, 1}^ where (U) 



1 T]{j) = 

otherwise. 



For a function / in Lp'{mN,Ki,...Kr) satisfying {f) N,Ki,...Kr = Oi define /o as the conditional 
expectation of / with respect to the a-field generated by vro. Then, similar to the proof 
of Theorem I3■14^ we project / on this u-field and on its orthogonal: 

(3.17) (/^) = ((/-/or) 

We first consider the second term. The arguments are similar to the ones used in the 
second step of the proof of Theorem 13.141 

We may think /o as a function defined on {0, 1}^. The generator acting on {0, 1}^ 
is the generator of the usual uniform symmetric exclusion process for which Quastel [8] 
proved a spectral gap. Therefore, we have 

We now turn to the first term of (|3.17|) . For a subset B C Ajv := {1, •••,-^} such that 
= Kq^ define /o,_b as a function on as follows: 



In this formula, stands for the configuration r] G such that 

Jo if i G B 

+ i if m=j- 

With this notation, we rewrite the first term of ()3.17p as follows: 

((/-/or) 

= ^ mN,K^,...Kr{HTl{i) = for all i G -B})(/o,B)mjv-Ko,Ki,...K,_r 
BcAjv 

The generator acting on ^^j^j^^b Ki K2 Kr-i ^® operator L^pf^Xo ■ 

induction assumption, there exists some constant Cr~i such that 



T]{i) 



(/o,B)mjv_Ko,ifi,...if,_i ^ Cr-l{L]^lKohB,fo,B 



i^n-Kq,Ki, 



N 



By the assumption Kq < therefore we have < and 

(/o,B)mjv_Ko,Ki,...K^_i — ^r{Liyf,f)mN,K^,...Kr 

for some constant Cr- □ 
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3.5 Closed Forms 

In this section, we introduce tlie notion of closed forms associated to our model and those 
of the generalized exclusion process. We have one-to-one correspondence between them, 
so algebraic characterization of the closed forms can be reduced to that for the generalized 
exclusion process. 

Let 'Hx,xJrei be a subspace of x = {—1,0, l}^*^ such that 
W,,,+e. := {r? ; + + *-rr(^) + n'-f'l^) + = !}• 

For two configurations 77 and ^ G Xi define D{ri,^) as follows: D(rj,^) = 1 if there 
exists a unique point x £ Tj'^ and a unique direction i such that r]{z) = ^(z) for all 
z ^ x,x + ei and r]{x) + ry(x + Cj) = ^{x) + ^{x + Cj) and r] ^, and D{r], £,) = otherwise. 
A path r(77, = iv = V^i ■•■) rf^''^,rf" = ^) from 77 to ^ is a sequence of configurations 
rf such that every two consecutive configurations satisfies D{rf' ,rf~^^) = 1. 

Consider a family of continuous functions u = (u^)i<i<d ^igz'* where : Tix^x+et ~^ ^■ 
For an ordered pair (?7,^) satisfying D{rj,$,) = 1, define a path integral -^(^^) of u by 
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A path integral can be naturally extended to paths of any length as 

m— 1 
j=0 

A family of continuous functions (u^)i<i<d,a;gz<i is called an /^-closed form if for all 
closed path r(r7,^), -^r(r)^)(") ~ ^ where a path r(r7,^) is called closed if 77 = ^. 

Next, let us recall the path integral and the closed form of the generalized exclusion 
processes. Let 'Hx,x+ei be a subspace of x = {—1, 0, 1}^'* such that 

n^,^ := {77 ; ^tf,+^^ir^) + ^^o'Tiri) + *f-|'''(77) + *S;o+^'(77) = 1}. 

For two configurations r] and ^ G Xi define D{r],(,) as follows: D{rj,$,) = 1 if there 
exists a unique point x G and a unique direction i such that ?7(z) = ^(z) for all 
z 7^ x,x -|- ej and 77(x) — 1 = S,{x) and 77(x + Ci) + 1 = ^{x + e^) or ?7(x) + 1 = ^(x) 
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and r]{x + e^) — 1 = ^(a; + e^), and D{r],^) = otherwise. A path T{r],$,) = {rj = 
77^, 77^, rf^-'^^rf^ = ^) from 77 to ^ is a sequence of configurations rf such that every two 
consecutive configurations satisfies D{rf,if'^^) = 1. 

Consider a family of continuous functions {ul;)i<i<d,xez<'' where : Hx,x+ei For 
an ordered pair (77, ^) satisfying D{ri,^) = 1, define a path integral if^^^ by 

j2 (~\ ._ i^iiv) if - 1 = .^(x) and + gj) + 1 = ^(x + gj) 
{v,0^ ■ |-u].(^) if ?7(x) + 1 = ^(x) and ??(x + e^) - 1 = ^(x + Ci). 

A path integral can be naturally extended to paths of any length as 

m—l 

A family of continuous functions (u^)i<i<(i,a;ezd called an /^-closed form if /|^^^^(u) = 

hold for all closed paths r(?7,^). 

Now, we construct one-to-one map from the set of /^-closed forms to the set of 7^- 
closed forms. For an /^-closed form (u^)i<j<rf ^.g^d, define a family of continuous functions 



4(77) := -^<(r?)*?,r^'(r7) + -^<(77)*S;!r^ 



{ri) 



then, (u^)i<i<<i,a:ez<* is an /^-closed form. On the other hand, for an /^-closed form 
{K)i<i<d,x&d., define a family of continuous functions {^l)i<i<d,x&d as 



+ v^ol(7,)^^;!r^(r7) + ^/Cc«l(r?)M^S;o+^^(r,) 



^~i^^^x=l,x+ei=-l^^ + b1(77^=°'^+"*=°))] *^fj''(r?), 

then, {'o'x)i<i<d,x&d is an /^-closed form. 

Let us introduce the notion of an /^-germ of closed form and an 7^-germ of closed 
form. A family of continuous functions {Qi)i<i<d (resp. Qi) where Qi : Ho^ei — M is an 
/^(resp. /^)-germ of closed form if := TxQi is an /^(resp. J^)-closed form. 

Main theorem of this section is formulated as follows: 

Theorem 3.18. For every I^-germ of closed form Q, there exists a sequence of L?{yp)- 
functions hn and constants (ci)i<j<d such that 

d 

Qi = Irni (cj '^{W)i + Vo^ei^hn) ^'^ L^i^p) for all 1 < i < d. 
j=i 
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By the one-to-one correspondence between the I^-germ of closed form and the I^- 
germ of closed form and Cauchy-Schwartz inequality, in order to prove Theorem 13.181 we 
have only to prove the next theorem: 

Theorem 3.19. For every I'^-germ of closed form g, there exists a sequence of Lp'{vp)- 
functions hn and constants (ci)i<j<rf such that 

d 

Qi = lim {a y^(l(^)i + Vo,eThJ in L?'{yp) for all 1 < i < d. 
i=i 

where Ui{r]) = l{,,(o)>o,r,(e,)<o} "■'^d 

Vo,eJO?) = <^'(r?)(/(r?°'^0 - /(^)) + ^l%{v){f{ri'''^) - /(r/)) 

+ <l\(^)(/(^°=°'^'=°) - fiv)) + Kfiiv){f{v'=-''''=') - fiv))- 

Applying the method in [3] , we can deduce this theorem since the proof depends only 
on the spectral gap estimates, which is proved by Theorem 1.8.1, and the fact that fp is 
translation invariant. 



4 Proof of Theorem 12.2 



The strategy of the proof is the same as given for liquid-solid system in [2j . The main step 
is to establish the local ergodic theorem. 

First, we consider a class of martingales associated with the empirical measure as 
similar to Case 1. We also use the same notations: T, M^{t), vr^, Q^n and Q* which are 
defined in Section 3. In Case 2, we assume the gradient condition, C++C- — Ca — '^Ce = 0, 
so the martingale M^{t) is rewritten as 

M^{t) = (vrf - {7r^,H) - f N''^ i^'' H){^)rMVs)ds 

where /i is a cylinder function defined by 

Hv) = C'+l{^(o)=i} - C_l{^(o)=_i} 
and A'^ H represents the discrete Laplacian of H: 



1=1 

Following the same argument in Section 3, it is easy to prove that 



lim P iv[ sup \M^{t)\ >6]=0 



N^oo ^ 'o<t<T 



for every 6 > 0, the sequence {Q^jv,A^ > 1} is weakly relatively compact and that every 
limit points Q* is concentrated on absolutely continuous paths iTt{du) = ii{t,u)du with 
density bounded by 1 and -1 from above and below respectively: — 1 < 7r(t, u) < 1. 
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From Theorem 16.11 stated below, there exists at most one weak solution of (12. 5p . 
Therefore, to conclude the proof of the theorem, it remains to show that all limit points 
of the sequence {Q^n,N > 1} are concentrated on absolutely continuous trajectories 
■K{t,du) = TT(t,u)du whose densities are weak solutions of the equation (|2.5p . For this 
purpose, all we have to show is that 

hm E^^[ TiV-'^ Yl i^''H)QrMris)ds] = Eq.[ f [ (A/7)(tz)P(7r(s, 
Jo ^ ^ Jo J-f 



where P is the function defined in ()2.6p . 

First, we establish the local ergodic theorem which enables us to replace the sample 
mean of microscopic variables with their average under the equilibrium measure having a 
microscopically defined sample density as its density-parameter. Let /i^ G ^(Xat) be the 
probability distribution of rjt on Xn /^^ be the space-time average of {^t^}o<t<T 

defined by 

T 



Then, we have that 

Proposition 4.1 (local ergodic theorem). For every cylinder function f, 



lim limsupS^'^[|/o,i^(r/) - (/)^x(o)|] = 
where Jq^k = (2K+ir T.\x\<K^xfir]). 

Proof. Let L be an operator on C defined by {Lf){ri) = '}2b&{'L<i)* ^bfiif)-, where {U^)* 
stands for the set of all directed bonds. Following the method of the proof of Theorem 4.1 
in [2], it is easy to show that {^l^^n is tight in V{x'^) and an arbitrary limit ^ G P{x^) 
satisfies f^{Lf) = for every cylinder function /. Moreover, by definition, fj, is invariant 
under spatial translations. Therefore, we have that support (/i) C {—1, 0}^'* U {0, 1}^'' . It 
is known that the translation-invariant L-stationary measure on { — 1,0}^'* or {0,1}^'* is 
a superposition of Bernoulli product measures. Then, the law of large numbers concludes 
the proposition. □ 

Next, we need to prove that the sample density defined microscopically can be re- 
placed in the limit with the macroscopic one. We can use Young measures to complete it 
by the exactly same way as in [2]. Then, combining with these results, the main theorem 
will be concluded. The details are omitted. 



5 Proof of Theorem 12.3 

In Case 3, since the hydrodynamic equation is the heat equation, the replacements which 
are required in Case 2 are unnecessary. So, this is the easiest case to prove the hydrody- 
namic limit. 
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First, we consider a class of martingales associated with empirical measure again. 
Then, as in Case 2, the martingale (t) is rewritten as 

M^(t) = {7Tf,H) - {7r^,H) - Ce [\7T^,A''H)ds 

Jo 

where A^ff represents the discrete Laplacian of H. 

Following the same argument in Section 3 and 4, all we have to show is that 

(5.1) Jim E^n[ ['n-^ Y1 hvs(x)=o}ds] = 0. 

With the notation defined in the last section, (jS.ip is rewritten as 

hm ii;^'~[l|,(o)=o}] = 0. 

Following the method of the proof of Theorem 4.1 in [2] again, it is easy to show that 
{p^}n is tight in V{x'^) and an arbitrary limit p S V{x'^) satisfies fJ^iLf) = for every 
cylinder function /. Moreover, by definition, /i is invariant under spatial translations. 
Therefore, we have that support (p) C {—1, 1}^ and it concludes the theorem. 



6 Weak Solutions of Nonlinear Parabolic Equations 

In this section, we fix the terminology of weak solutions of parabolic equations and present 
the uniqueness of such equations. Hereafter : M — > M is a strictly increasing and Lipschitz 
continuous function. We consider the Cauchy problem: 

. dtp{t, u) = A{cl>{p{t, u))) ( = Q^^^Pi^^ ^))) 

/5(0,-) =Po(-), 
and define weak solutions of this Cauchy problem. 

Definition 6.1. Fix a bounded initial profile /jq : TT'^ — > M. A measurable function p = 
p{t) = p{t,u) G C{[0,T],M{T'^))nL'^{[0,T]xT'^) is a weak solution of the Cauchy problem 
116. 1\) if for every function H -.T'^ of class C^(T'^) and for every < t <T 

H{u)p{t,u)du = j ds (t){p{s,u))AH{u)du+ / H{u)po{u)du. 



/Jd. Jo JT'' Jjd. 

We prove the uniqueness of weak solutions in this class. 

Theorem 6.1. Fix a bounded measurable function po ■.T'^ ^ M. There exists at most one 
weak solution of the parabolic equation 16. 

Proof. We can apply the proof of Theorem A. 2. 4. 4 in ^ straightforwardly. □ 
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